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ABSTRACT 

We present an optical multicolor-imaging study of the galaxy cluster Abell 
521 at z — 0.25, using Suprime-Cam on the Subaru Telescope, covering an area 
of 32 x 20 arcmin 2 (9.4 x 5.8h^ Mpc 2 at z = 0.25). Our imaging data taken 
with both a narrow-band filter, NB816 (A = 8150A and AA = 120A), and 
broad-band filters, B, V,Rc,i', and z' allow us to find 165 Ha emitters. We 
obtain the Ha luminosity function (LF) for the cluster galaxies within 2 Mpc; 
the Schechter parameters are a = -0.75 ± 0.23, (p* = io _0 - 25±0 - 20 Mpc" 3 , and 
L* = io 42 03±017 erg s -1 . Although the faint end slope, a, is consistent with that 
of the local cluster Ha LFs, the characteristic luminosity, L*, is about 6 times 
(or « 2 mag) brighter. This strong evolution implies that Abell 521 contains 
more active star-forming galaxies than the local clusters, being consistent with 
the observed Butcher-Oemler effect. However, the bright L* of Abell 521 may 
be, at least in part, due to the dynamical condition of this cluster. 

Subject headings: galaxies: clusters: individual (Abell 521) — galaxies: evolution 
— galaxies: luminosity function, mass function 
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1. INTRODUCTION 

Recent observations of clusters of galaxies at intermediate-redshift (0.1 < z < 0.5) 
demonstrated that environmental processes play an important role on both the star-formation 
activity and the morphology of galaxies (e.g., Dressier et al. 1997; Balogh et al. 1998; Pog- 
gianti et al. 1999; Moss & Whittle 2000). The galaxy evolution in such clusters may be 
explained as the result of the increased activity of the field galaxies at more distant redshift 
(Lilly et al. 1996), modulated by changing the infall rate onto the clusters (Bower 1991). 
However, little information on the global star formation history in such intermediate-redshift 
clusters has been obtained even to date. 

In many previous spectroscopic studies such as the Canadian Network for Observational 
Cosmology (CNOC) and MORPHS surveys (e.g., Abraham et al. 1996; Poggianti et al. 
1999), the [O ll]A3727 emission line has often been used to investigate the star formation 
activity in cluster galaxies. However, the Ha emission line provides a more reliable indicator 
of star formation because it is less sensitive both to dust extinction and to metallicty effects 
(Kennicutt 1998). Recently, an Ha survey has been made for nearby clusters of galaxies 
Abell 1367 and Coma (Iglesias-Paramo et al. 2002). These data are utilized to derive 
the Ha luminosity function (LF) of the cluster member galaxies. Optical spectroscopic 
surveys are also useful in studying Ha emission-line properties of cluster member galaxies 
(Balogh et al. 2002; Couch et al. 2001). However, such spectroscopic sample are selected 
by using optical broad-band color properties, being independent from their Ha emission-line 
properties. Therefore, these data sets may miss a part of Ha-emitting galaxies that are too 
faint to be selected in a continuum-magnitude-limited sample. In order to determine the 
shape of the Ha LF unambiguously, it is required to carry out a deeper Ha imaging survey. 

Motivated by this, we performed our deep Ha imaging survey of a cluster of galaxies, 
Abell 521 at z = 0.25 which is a relatively rich (richness class 1) cluster; note that its Bautz- 
Morgan Type is III (Abell, Corwin, & Olowin 1989). Recent studies of this cluster showed 
that the projected galaxy density distribution in the central region has a very anisotropic 
morphology, as it exhibits two high-density filaments crossing in an X-shaped structure at 
the barycenter of the cluster (Arnaud et al. 2000; Maurogordato et al. 2000; Ferrari et 
al. 2003). The observed high velocity dispersion of the cluster (w 1325 km s -1 ) may be 
attributed to the presence of those substructures. Since this velocity dispersion is also high 
compared to that expected from the temperature of the X-ray gas (kT = 6.3 keV, Arnaud 
et al. 2000), it is suggested that this cluster is undergoing strong dynamical evolution driven 
by mergers among substructures. In this paper, we present observational properties of Ha- 
emitting galaxies in Abell 521 and then derive the Ha LF for the first time. Throughout this 
paper, magnitudes are given in the AB system (Fukugita, Shimasaku, & Ichikawa 1995). We 
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adopt H = 50 km s 1 Mpc 1 , q = 0.5, and A = cosmology, for comparison with previous 
results in the literature. 



2. OBSERVATIONS AND DATA REDUCTION 

The imaging data of Abell 521 were obtained with Suprime-Cam (Miyazaki et al. 2002) 
on the 8.2 m Subaru Telescope atop Mauna Kea during two runs (2001 October and 2002 
February). We used the narrow-band filter, NB81Q (A = 8150A and AA = 120A); the 
central wavelength corresponds to a redshift of ~ 0.24 for Ha emission. We also used broad- 
band filters, B, V, Rc, i', and z'. A journal of our observations is summarized in Table 
1. 

We used IRAF and the mosaic-CCD data reduction software, NEKOSOFT (Yagi et al. 
2002), to reduce and combine the individual CCD data. Photometric and spectrophotometric 
standard stars used for the flux calibration are as follows; SA98 (Landolt 1992) for the B 
data, SA95 and SA98 (Landolt 1992) for the V data, SA95, SA113, and PG0231 (Landolt 
1992) for the R c data, GD50 (Oke 1990) for the i' data, and GD50, GD108, and PG1034+001 
(Stone 1996) for the NB816 data. Since the zero point of Landolt (1992) is based on Vega, 
we convert the B, V, and Rc magnitudes to the AB system assuming monotonic colors of 
5(Vega) - B(AB) = 0.140, V(Vega) - V(AB) = 0.019, and i? c (Vega) - R C (AB) = -0.169 
(Fukugita et al. 1995). For i' and NB816, we multiplied the atmospheric, instrumental, 
and filter transmittance, mirror reflectivity and CCD quantum efficiency with the SED 
of standard stars to calculated absolute flux. The z' data were calibrated by using the 
magnitude of the quasar SDSSp J104433.04-012502.2 (Fan et al. 2000), which were obtained 
in the same observing night (see Fujita et al. 2003). Since any quasar is a potentially variable 
object, the photometric calibration of the z' data may be more unreliable than those of the 
other band data. In order to estimate the uncertainty, we compared the colors of stars in 
our area with Galactic stars given by Gunn & Stryker 7 and applied the correction to our 
z' data. The combined images for the individual bands were aligned and smoothed with 
Gaussian kernels to match the PSF (= I'M). The final images cover an area of 32 x 20 
arcmin 2 , which corresponds to 9.4 x 5.8 Mpc 2 at z — 0.25. 

Source detection and photometry were performed by using SExtractor version 2.2.2 
(Bertin & Arnouts 1996). The NB&IQ image was used as a reference image for detection, 
and photometry was performed for each band image. We used a criterion that a source has at 
least a 13-pixel connection at 2a level. For each object, we measured the magnitude within a 



7 Bruzual-Persson-Gunn-Stryker (BPGS) Library (ftp://ftp.stsci.edu/cdbs/cdbs2/grid/bpgs) 
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fixed aperture of a 3" 15 kpc at z = 0.25) diameter for colors and the magnitudes in all 6 
band images. However, we used an adaptive aperture of diameter 2.5rx for measurement of 
total magnitudes, where r K is the Kron radius (Kron 1980). Our 3cx rms within a 3"0 aperture 
in each band is given in Table 1. 

In order to correct for the Galactic extinction, we calculated the extinction at each 
band (Cardelli et al. 1989), assuming that A B = 0.324 (Schlegel et al. 1998, NED 8 ). The 
extinction at each band, A\, is given in Table 1. 



3. SELECTION OF Ha EMITTERS 

We defined the off-band continuum flux of objects as i'z' = 0.6i' + OAz', which was 
computed from a linear interpolation between the effective frequencies of the i' and z' bands. 
Taking the scatter in the i'z' — NB816 color into account, emission-line objects were selected 
with the following criteria, 

i'z'-NB8l6 > max (0.1, 3(7 error of i'z'- NBS16), (1) 
i'z' < 24.9. (2) 

The color of i'z' — NB816 = 0.1 corresponds to an equivalent width of 12A in the rest frame 
at z — 0.247. The magnitude limit corresponds to the 3<j rms of the i'z' band. These criteria 
are shown by the solid, dashed, and dotted lines in Figure 1. Then we found 324 objects 
that satisfy the above criteria; note that they are brighter than the limiting magnitude at 
each band. 

A narrow band imaging can potentially detect galaxies with different emission lines 
at different redshifts. Emission lines detected in our narrow-band imaging survey are, for 
example, Ha (z « 0.24), [O ill] AA4959, 5007 (z « 0.63), H/3 (z « 0.68), and [O ll]A3727 
(z R3 1.18). In order to distinguish Ha emitters at z ~ 0.24 from other emission-line objects, 
we investigate their broad-band color properties. In Figure 2, we show the diagram between 
B — V and V — i! for the 324 emission-line objects. We also show the model galaxy colors 
that estimated by using the population synthesis model GISSEL96 (Bruzual & Chariot 1993). 
The star formation rate of model galaxies is proportional to exp(— t/r), where we use r = 1 
Gyr. Changing the age of model galaxies, we generate the SEDs of five types of galaxies [t = 
1, 2, 3, 4, and 8 Gyr for SB (young starburst), Irr, Scd, Sbc, and E, respectively]. Then we 



8 NASA/IPAC Extragalactic Database http://nedwww.ipac.caltech.edu/ 
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find that Ha emitters at z ps 0.24 can be selected by adopting the following color criterion, 

B — V > 0.75(V — i') + 0.25. (3) 

As shown in Figure 2, this criterion appears to isolate Ha emitters from other emitters at 
different redshifts. We then obtain a sample of 165 Ha emitters. 

In order to examine whether or not this color criterion works well more quantitatively, we 
investigate color properties of the spectroscopically confirmed member galaxies in Abell 521. 
The previous spectroscopic surveys identified 125 member galaxies in this cluster (Ferrari 
et al. 2003; Maurogordato et al. 2000). These members are identified in our catalog. We 
found that 123 objects satisfy the color criterion based on our photometry. The remaining 
two objects were rejected in our selection procedure because they have neighboring galaxies, 
causing large photometric errors. As a further check, we examined the spatial distribution 
of the galaxies which are selected by the above color criterion. In Figure 3, the excess of 
surface number density by the cluster members is clearly seen in the area within the radial 
distance of 2 Mpc from the center of Abell 521; note that the center is taken at the position 
of the main X-ray cluster (Arnaud et al. 2000). 

In order to investigate the detection completeness of Ha emitters, we performed a detec- 
tion test on simulated artificial images. We prepared model galaxies with broad-band colors 
of late type galaxies (SB, Irr, Scd, and Sbc) at z — 0.24 generated by GISSEL96 model. 
Since the equivalent width of Ha emission is generally different form galaxy to galaxy, we 
put the excess of NB816 to i'z' continuum as a parameter. We studied a color range of 
0.1 < i'z'-NB81Q < 1.1 with five bins of A (i'z 1 - NB81Q) = 0.2. Since the total magnitude 
of A^_B816 is another parameter, we studied a magnitude range of 18.0 < NB81Q < 24.0 
with twelve bins of A(iV.B816) = 0.5. We prepared 500 model galaxies in each set of NB 816 
and i'z' — NB816. The light distribution of model galaxies is generated as the exponential 
law by IRAF ARTDATA. Their sky positions, half-light radius (1 to 7 kpc), and ellipticities 
(0.3 to 1.0) were randomly set and put into the CCD data together with Poisson noises. 
After smoothing model-galaxy images to match to the seeing size, we tried to detect them 
using SExtractor and to select as Ha emitters with the same procedure as that used for the 
observed data. The detectability of the model emitters is shown in Figure 4 as a function 
of NB816 magnitude. We found that the detectability is higher than 80% for objects with 
NB816 < 22.0, but lower than 30% for objects with 23.0 < iV£816 < 24.0. 

Next, we examined the contamination by emission-line objects at other redshifts (e.g., 
[O ill] AA4959, 5007 at z ps 0.63) in the same manner as that for the detection completeness 
estimate. We found the contamination to be only less than 5% even for the faint objects 
with 23.0 < NB816 < 24.0. In conclusion, our color criterion allows us to select reliable Ha 
emitters at z 0.24 for sources with NB81Q < 23.0. 
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4. Ha LUMINOSITY FUNCTION 

4.1. Emission-Line Equivalent Width 

At first we derive the emission-line equivalent width of the individual Ha emitters. Since 
the pass band of our narrow-band filter is too wide to separate [N ll]AA6548, 6584 from Ha, 
the equivalent with derived here is that contributed from the three emission lines. We adopt 
the same method as that used by Pascual et al. (2001) to calculate Ha + [N ll]AA6548, 6584 
(hereafter Ha + [N n]) equivalent width. The flux density in each filter can be expressed as 
the sum of the line flux and the continuum density (the line is covered by both filters), 

Inb = fc + (F L /ANB), (4) 
/iv = /c + 0.6(F L /Ai'), (5) 

where fc is the continuum flux density, the line flux, ANB and Ai' the NB816 and i' 
band filter effective widths, respectively, and J'nb and fi> z > the flux density in each filter. 
Then the equivalent width can be expressed as follows; 



EW ohs (Ra + [N II]) = ^ = ANB 

Jc 



Inb — fi' 



Ji>z> - 0.6f NB (ANB/Ai>) 
where Ai' = 1535A for the i' band and ANB = 120A for the NB81Q band. 



(6) 



4.2. Ha Luminosity 

In order to obtain the Ha luminosity, we have to correct for the flux contribution from 
the [N n] doublet. We also have to apply correction for internal extinction. For the two 
corrections, we adopt the flux ratio of Ha/[N il] = 2.3 (Kennicutt 1992) and A Hq = 1 mag as 
a typical internal extinction at z ~ 0.2 (Tresse & Maddox 1998). Applying these corrections, 
we estimate the Ha luminosity, L(Ha), from the total magnitude for each objects. We then 
obtain the star formation rate, SFR, using the following relation (Kennicutt 1998), 

SFR (M© yr" 1 ) = 7.9 x 10" 42 L(Ha) (erg s" 1 ). (7) 

The net Ha flux depends on the source redshift because of the fixed transmittance of 
our NB816 filter. In order to correct for this effect, we adopt the following assumptions. (1) 
The velocity (redshift) distribution of the cluster galaxies is a Gaussian probability density 
function P(z) around a mean velocity of 74019 km s~ x with a dispersion of 1325 km s^ 1 
(Ferrari et al. 2003). (2) The equivalent width distribution of the Ha emitters in the rest 
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frame is expressed as the following exponential functional form; Q(EW) = Aexp(— EW/B), 
where A and B are free parameters. Then, the rest frame equivalent width probability of 
the Ha emitters at redshift z can be expressed as P{z) x Q(EW). The best-fit parameters 
are obtained so as to fit the distribution of i'z' — NB81Q. 

We apply the correction for the incomplete velocity coverage of the A^_B816 filter. Since 
the number of the Ha emitters excluded from the filter can be calculated from the value 
of the equivalent width, we correct the equivalent width distribution function that derived 
above. Therefore, using the equivalent width distribution function, we obtain the corrected 
Ha luminosity distribution. 

4.3. Statistical Weights 

Since the excess of surface galaxy density by the cluster member is clearly seen in the 
area within 2 Mpc radius (see Figure 3), it seems safe to use the Ha emitters located in 
the region; note that there are 53 Ha emitters in this region. Although the NB816 filter 
corresponds to the redshifted Ha emission line at z — 0.242 ± 0.009, the mean redshift of 
the cluster is z = 0.247. Therefore, our Ha emitter sample may suffer contamination from 
field Ha emitters, in particular from field galaxies in the foreground. In order to correct for 
this effect, we estimate the field contamination by the following method. 

In order to obtain a sample of field galaxies, we use an area of 3 Mpc outside from the 
cluster center, taking account that the virial radius is estimated to be 2.74 Mpc (Girardi & 
Mezzetti 2001). However, since a subclump is seen to the North-East at about 4.5 Mpc from 
the cluster center, we exclude this area from our analysis. The remaining sky area is used 
to estimate the surface number density of field galaxies as a function of A^_B816 magnitude. 
Since Abell 521 shows elongated galaxy distribution (Ferrari et al. 2003), we use the local 
surface number density, which is calculated from 10 nearest neighbors, rather than the radial 
surface number density distribution. In this way, all Ha emitters in the cluster region are 
weighted by the luminosity-dependent factor, w\ and the density-dependent factor, u>d, after 
renormalizing Wd so that the sum of the total galaxy weights, W\ x u? d is equal to the sum of 

W\. 

4. 4. Ha Luminosity Function 

We derive the Ha luminosity function (LF) for the cluster region. In order to normalize 
the Ha LF to a proper volume, we use the sphere volume within 2 Mpc radius. The Ha LF 
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for our sample is well fitted by the Schechter (1976) function; 

0(L)dL = (j)\L/L k ) a exp(-L/L*)d(L/L*). (8) 

The size of the bins is taken to be AlogL(Ha) = 0.5. The number counts in the individual 
bins are given in Table 2, and Figure 5 shows the Ha LF for Abell 521 and its Schechter 
fit. After correcting for the detection completeness (see section 3), we obtain the best-fit 
parameters for the cluster membership within 2 Mpc radius are 

a = -0.75 ±0.23, 
0* = 10-°- 25±a2 ° Mpc" 3 , 
L* = 10 42 - 03±a17 erg s" 1 . 

The inset panel in Figure 5 shows the \ 2 error contours corresponding to the 1 and 2a levels. 
It can be seen that a and L* values are highly correlated. 

5. DISCUSSION 

In Figure 5, we also show the Ha LFs of field galaxies in the local Universe (Gallego et 
al. 1995), at z ~ 0.2 (Tresse & Maddox 1998), at z w 0.24 (Fujita et al. 2003), and those 
of the nearby clusters of galaxies (Iglesias-Paramo et al 2002). As for the nearby clusters 
of galaxies, Iglesias-Paramo et al. (2002) derived the Ha LFs for Abell 1367 and Coma 
(richness class 2). The physical areas analysed by them are nearly the same as that of Abell 
521. Although they applied a slightly different extinction correction; A Ha = 1.1 mag for 
type Scd or earlier and A Ua = 0.6 for type Sd or later (Boselli et al. 2001), they obtained 
a ~ —0.7 and L* pa 10 41 ' 25 erg s _1 for these clusters. The inset panel in this figure shows 
the best-fit parameters of Abell 521 and the others. There is the significant difference in 
Ha LFs between the clusters and the fields. It is found that the faint end slope, a, of the 
Ha LF of Abell 521 is flatter than that of the field Ha LFs at the same redshift. Therefore, 
although the sample analysed here is not so large, it appears that the flatter faint end slope 
of the Ha LF is commonly seen in the clusters of galaxies between z ~ and z ~ 0.2. 

The most important finding in this study is that the characteristic luminosity, L* of 
Abell 521 is ~ 6 times (or pa 2 mag) brighter than those of the local cluster Ha LFs, and 
nearly the same as those of the field Ha LFs. The evolution of the field Ha LF from z ~ 0.2 
to z ~ appears closely related to the 0* parameter rather than to the L* and a parameters. 
However, as for the cluster Ha LFs, the L* parameter appears to evolve strongly from z ~ 0.2 
to z ~ 0. Furthermore, we investigate the following galaxy clusters for which spectroscopic 
Ha emitter surveys were already made; Abell 1689 at z — 0.18 (Balogh et al. 2002) and AC 
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114 at z — 0.31 (Couch et al. 2001). Since their observations are based on spectroscopy, 
it is difficult to compare with our Ha LF straightforwardly. Therefore, we compare the Ha 
LF of Abell 521 with the other clusters, only using the Ha luminous galaxies (L* > 10 415 
erg s _1 ). It should be noted that the same dust extinction correction (i.e., An a = 1 mag) is 
applied for these clusters, and the Ha luminosity from spectroscopy may be underestimated 
due to aperture bias. 

Figure 6 shows the surface density of Ha luminosity distribution of each cluster; note 
that the surface density is used instead of the volume density which is used in Figure 5, to 
compare simply with previous results. Although there are some galaxies whose SFRs exceed 
SFR ~ 4M yr _1 in Abell 521 and Abell 1689, no such galaxies are found in Abell 1367 
and Coma. This SFR is about the same level as that observed in the Milky Way (Rana 
1991). This result implies that Abell 521 and Abell 1689 contain more active star-forming 
galaxies than the nearby clusters. However, AC 114 has nearly the same property in the Ha 
LF as those of the nearby clusters in spite of its higher redshift. 

This evolution may simply reflect the Butcher-Oemler (BO) effect, which is the increase 
in the fraction of blue galaxies, presumably star-forming galaxies, in clusters with increasing 
redshift (Butcher & Oemler 1984). In fact, the fraction of blue galaxies in our Ha emitters 
by the definition of Butcher & Oemler (1984) is ~ 85%. The blue fraction, /b, in Abell 
1689 is comparable to that of Abell 521 (/b ~ 0.17, see for details Appendix) and much 
higher than those in the nearby rich clusters such as Coma (/b = 0.03). However, AC 114 
and Abell 1367 also have a large fraction of blue galaxies (/b ~ 0.20). Therefore, the star 
formation activity in clusters of galaxies may not be understood solely by the BO effect. 

These results suggest that some physical processes work in their star formation activity, 
independent of the blue fraction of each cluster. We investigated what parameter of cluster 
correlates with the existence of high SFR galaxies. One possible parameter is the dynamical 
status, as pointed out by Balogh et al. (2002). Although Abell 1689 appears to show a round 
shape in X-ray, the observed radial velocity distributions of the cluster shows substructures. 
On the other hand, AC 114 appears to be relaxed dynamically, although the cluster member 
galaxies shows an elongated morphology (Balogh et al. 2002). Therefore, it seems important 
to investigate the radial velocity distributions for the Ha emitters of each cluster (Figure 
7). In this figures, we use the relative radial velociy, Av r , to the systemic velocity, t> r ,o, of 
each cluster, and t> r ,o is shown in the left corners of each figure. For Abell 521, Abell 1367, 
and Coma, their redshift are unknown from their imaging survey, since the Ha emitters 
are detected using the narrow-band filters. Therefore, their distributions are shown using 
available spectroscopic data in the literature (Ferrari et al. 2003; Iglesias-Paramo et al. 2002; 
Cortese et al. 2003). The vertical dashed lines in this figures show the velocity coverage 
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corresponding to the FWHM of the narrow-band filters. For Abell 521, in particular, since 
some Ha emitters with At> r > km s _1 may be excluded from the filter transmittance, we 
also show the velocity distribution of galaxies identified with other emission lines ([O II] or 
H/3) taken from Ferrari et al. (2003). It is found that this radial velocity distribution is 
consistent with that of the Ha emitters with Av r < km s _1 . Therefore, we expect that it 
also traces that of the Ha emitters with Av T > km s _1 . In this figure, the radial velocity 
distribution for the Ha emitters in Abell 521 is inconsistent with a Gaussian distribution 
which we assumed in section 4.1. However, since the spectroscopic sample is rather small, 
it is dangerous to apply this distribution to our sample. Even if this distribution is valid, it 
does not change the result that L* of Abell 521 is much brighter than those of the nearby 
clusters. 

In Figure 7, the radial velocity distributions in Abell 1367 and Coma tend to be concen- 
trated around Av r km s" 1 . Therefore, it is suggested that these galaxies are associated 
with a galaxy population in the cluster core, probably a virialized system. For AC 114, 
the radial velocity distribution suggests that the most Ha emitters are associated with the 
structure centered at At> r m +2000 km s -1 . Even if they comprise a substructure infalling 
on to the cluster center, the star-forming activity seems to have already been ceasing by the 
environment effect (Kodama et al. 2001; Gomez et al. 2003). Abell 521 and Abell 1689 show 
that the radial velocity distributions show several peaks in a wide velocity range, and there 
are few Ha emitters at the mean cluster redshifts. And for Abell 521, they are not thought 
to be field contamination but be associated with the cluster, since their number density is 
high at the center of the cluster (Figure 3). The Ha emitters in the two clusters are likely to 
be in the process of accretion of the field population. If this is the case, we may explain why 
their star-forming activities are higher than those of the other clusters. It is thus suggested 
that the star formation activity in clusters seems to be strongly related to the dynamical 
status of the clusters. 
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A. COMMENTS ON THE BLUE FRACTION 

In order to derive the fraction of blue galaxies, /b, in Abell 521, we followed the proce- 
dure of Butcher & Oemler (1984). We calculated the fraction of galaxies that are 0.2 mag 
bluer than the color-magnitude sequence of early-type galaxies in the rest-frame B — V color 
within a magnitude limit, My = —20, and that are enclosed within R 30 , the radius which 
contains 30% of the cluster population. 

The blue fraction, / B , was calculated by subtracting the number count expected for the 
field from that obtained for the cluster sample. In order to obtain the value of / B , we used 
the V — Rq color index, which is very close to B — V rest frame at the cluster redshift. The 
magnitude limit, My = —20, in the rest frame corresponds to Rq — 20.8 at cluster redshift, 
and a difference of A(B — V) = 0.2 mag is A(V — R c ) — 0.17 mag. The value of i? 30 is 
determined from the surface number density profile (Figure 3). We found R^q = 3.1 arcmin 
0.9 Mpc), being very similar to the typical radius of R30 ~ 1 Mpc for a sample of rich 
clusters at 0.2 < z < 0.4 (Kodama & Bower 2001). As a result, we obtain f B = 0.17 ± 0.03. 
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Table 1. Journal of the observations 



Filter 


Obs. Date 


Total Exp. (s) 




Limiting Mag. b 


B 


2002 Feb. 15 


480 


0.324 


25.3 


V 


2001 Oct. 14, 15 


2250 


0.245 


25.5 


Rc 


2001 Oct. 14 


1620 


0.203 


25.5 


i' 


2001 Oct. 14 


1320 


0.157 


25.4 


NB81Q 


2002 Feb. 14 


3600 


0.139 


24.6 


z' 


2002 Feb. 15 


840 


0.113 


24.2 



a Galactic extinction are applied to our imaging data (Cardelli et al. 
1989; Schlegel et al. 1998). 

b The limiting AB magnitude (3a) within a 3"0 aperture. 



Table 2. Ha luminosity function 



logL(Ha) 


(erg s x ) 


Number of Ha emitters 


log0 (Mpc~ 3 ) 


40.5 - 


41.0 


10.62 


-0.60 


41.0 - 


41.5 


12.77 


-0.52 


41.5 - 


42.0 


11.81 


-0.55 


42.0 - 


42.5 


5.23 


-0.90 


42.5 - 


43.0 


0.17 


-2.40 
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i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 — :- | — — •• s — j 

!'z'-A«816 = 0.1 '. : V . 




NBS16 

Fig. 1. — The color-magnitude diagram of i'z' — NB816 vs. NB816. The horizontal solid 
line corresponds to i'z' — NB816 = 0.1. The dashed line shows the distribution of 3a error, 
and the dotted line shows the limiting i'z' magnitude. 
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V-V 

Fig. 2. — The two-color diagram of B — V vs. V — i' for all the emitters. The objects with 
i'z' < 20 are shown by open circles, those with 20 < i'z' < 22 are shown by filled triangles, 
and those with i'z' > 22 are shown by crosses. The error bars in the left corners of the 
plots indicate the mean uncertainties on colors. Colors of model galaxies (for SB, Irr, Scd, 
Sbc, and E) from z = (filled squares) to z = 2 (filled diamonds) are shown by blue lines. 
Colors of galaxies at z = 0.24, 0.64, 0.68, and 1.18 (for Ha, [O ill], H/3, and [O n] emitters, 
respectively) are shown with red, purple, yellow, and green lines, respectively. We select the 
objects above the dashed line as Ha emitters except for the two open squares. 
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Fig. 3. — Surface number density distribution of the galaxies with i'z' < 23 as a function 
of the radius from the center of Abell 521. The objects with B — V > 0.75(V — i') + 0.25 
are shown by the filled circles and those with B — V < 0.75(V — i') + 0.25 are shown by the 
open triangles. The vertical solid line shows the radius of 2 Mpc from the cluster center. 
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Fig. 4. — Detection completeness of the Ha emitters derived from the simulation. The 
horizontal line shows the detection completeness of 0.80. 
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Fig. 5. — The derived Ha luminosity function (LF) for Abell 521 is shown by the thick solid 
curve with the data points (open circles). For reference, the Ha LFs of Abell 1367 and Coma 
(Iglesias-Paramo et al. 2002) are shown by the thick dashed and dotted curves, respectively. 
We also show those of field galaxies in the local Universe (Gallego et al. 1995, thin dotted 
line), at z ~ 0.2 (Tresse & Maddox 1998, thin dashed line), and at z ~ 0.24 (Fujita et al. 
2003, thin solid line). Inset panel shows 1 and 2a error contours for the best-fit Ha LF 
parameters of Abell 521. The best-fit parameters of Abell 1367 and Coma are shown by the 
filled circle and the filled triangle, and those of field galaxies in the local Universe and at 
z ~ 0.2 are shown by the open circle and the open triangle, respectively. Note that those of 
field galaxies at z 0.24 are not shown in this figure. 
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Fig. 6. — The surface density of Ha luminosity distribution of the cluster galaxies for Abell 
521. We also show those for Abell 1367, Coma (Iglesias-Paramo et al. 2002), Abell 1689 
(Balogh et al. 2002), and AC 114 (Couch et al. 2001). 
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Fig. 7. — Radial velocity distributions for the Ha emitters of each cluster, with a binning of 
1000 km s~ x . We use the relative radial velociy to the systemic velocity, i> r ,o, of each cluster, 
and i> r ,o is shown in the left corners of each figure. The vertical dashed lines show the velocity 
coverage corresponding to the FWHM of the narrow-band filters for Abell 521, Abell 1367, 
and Coma. 



